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DEVELOPMENT  AND  TESTING  OF  THE  COUPLED 
OCEAN/ ATMOSPHERE  MESOSCALE 
PREDICTION  SYSTEM  (COAMPS) 

1.  Introduction 

U.  S.  Navy  operations  require  a  very  detailed  specification  of  the  environment. 
The  most  recent  example  of  this  was  exhibited  in  Operation  Desert  Storm.  The 
weapons  that  are  used  in  operations  such  as  this  are  affected  by  moisture,  temper¬ 
ature,  aerosols,  etc.  To  determine  these  parameters,  the  Navy  currently  relies  on 
the  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS)  (Hogan 
and  Rosmond,  1991)  and  the  Navy  Operational  Regional  Atmospheric  Prediction 
System  (NORAPS)  (Hodur,  1987).  Both  NOGAPS  and  NORAPS  are  hydrostatic 
atmospheric  models  using  horizontal  resolutions  in  the  range  of  40-150  km.  For 
operations  such  as  Desert  Storm,  the  Navy  needs  higher  temporal  and  spatial  reso¬ 
lution  than  either  NOGAPS  or  NORAPS  currently  provides.  Therefore,  the  Navy 
must  look  to  develop  models  capable  of  predicting  scales  of  motion  at  and  below 
the  meso-yS  scale,  which  will  require  grid  spacings  of  20  km  or  less. 

The  prediction  of  atmospheric  motions  of  meso-/?  and  smaller  scale  implies 
that  the  hydrostatic  approximation  may  be  invalid  at  times,  particularly  for  con¬ 
vection  and  smaller-scale  topographic  features  where  the  vertical  wavelength  is  a 
significant  fraction  of  the  horizontal  wavelength.  The  concept  of  using  the  nonhy¬ 
drostatic  form  of  the  equations  as  a  prediction  tool  started  with  the  work  of  Ogura 
and  Charney  (1962)  and  Ogura  and  Phillips  (1962).  However,  it  wasn’t  until  the 
1970’s  that  numerical  models  based  on  the  nonhydrostatic  equations  were  devel¬ 
oped.  Several  of  these  models  were  used  to  study  convective  processes  (Miller  and 
Pearce,  1974;  Schlesinger,  1975;  Clark,  1977;  Klemp  and  Wilhelmsor,  1978;  Clark, 
1979;  Tripoli  and  Cotton,  1982)  while  others  studied  mountain  flows  (Durran  and 
Klemp,  1983;  Schumann,  1987;  Xue  and  Thorpe,  1991).  The  British  Meteorologi¬ 
cal  Office  developed  a  nonhydrostatic  model  for  operational  use  (Tapp  and  White, 
1976;  Carpenter,  1979)  and  Tripoli  (1992a)  has  recently  developed  a  nonhydrostatic 


model  to  study  scale  interaction. 

One  of  the  problems  in  dealing  with  the  nonhydrostatic  equations  is  that 
they  allow  for  rapidly  propagating  sound  waves,  thereby  limiting  the  time  step  of 
the  model.  The  anelastic  approximation  (Ogura  and  Philhps,  1962)  eliminates  the 
sound  waves,  but  necessitates  the  solution  of  a  three-dimensional  elliptic  equation 
every  time  step  for  pressure  or  geopotential,  which  can  be  computationally  expen¬ 
sive  (Tapp  and  White,  1976;  Miller,  1974;  Clark,  1977;  Xue  and  Thorpe,  1991). 
Klemp  and  Wilhelmson  (1978)  used  the  fully-compressible  equations  with  a  time¬ 
splitting  method,  which  allowed  for  large  time-steps  to  be  taken  for  the  slow  modes 
and  shorter  time  steps  for  the  faster  sound  modes.  This  time-splitting  technique, 
coupled  with  semi-implicit  differencing  in  the  vertical,  made  the  computational 
expense  of  integrating  the  compressible  form  of  the  nonhydrostatic  equations  com¬ 
petitive  with  the  anelastic  form.  The  quasi- compressible  approximation  (Chorin, 
1967),  which  artificially  slows  down  the  sound  waves,  is  an  alternative  method  to 
allow  for  larger  time  steps  in  the  fully-compressible  equations  (Droegemeier  and 
Wilhelmson,  1987). 

The  recent  advances  made  in  mesoscale  modeling  coupled  with  the  products 
needed  by  the  Navy  suggest  a  redefinition  of  ”regioneil”  modeling.  It  is  now  apparent 
that  global  models,  currently  capable  of  resolutions  in  the  50-100  km  range,  can 
handle  synoptic- scale  systems  quite  well  (Hogan  and  Rosmond,  1991).  This  is 
the  range  in  which  regional  models,  such  as  NORAPS,  have  served  over  the  p2ist 
decade.  However,  due  to  the  history  of  armed  conflicts  in  the  vicinity  of  coastal 
areas  and  the  pace  of  computer  technology,  we  will  hereafter  refer  to  "regional”, 
or  "limited- area”  modeling,  as  the  development  of  a  data  zissimilation  system  for 
the  analysis  and  forecast  of  meso-/S  and  meso-7  scales  of  motion  in  coastal  areas. 
This  coastal  data  aissimilation  system  should  allow  for  nonhydrostatic  processes. 
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detailed  topographic  flows,  and  explicit  treatment  of  precipitation  processes,  and 
be  fully  coupled  with  the  lower  boundary,  with  detailed  description  and/or  coupled 
forecasts  of  surface  parauneters  such  as  the  sea-surface  temperature  (SST),  ground 
wetness  and  ground  temperature.  This  system  must  take  into  account  new  data 
sources  such  2is  Doppler  and  profiler  winds  as  well  as  newly  emerging  computer 
technology  such  as  parallel  computer  architecture. 

To  serve  the  Navy’s  future  mesoscale  modeling  needs,  the  Naval  Research  Lab¬ 
oratory  (NRL)  is  developing  the  Coupled  Ocean/ Atmosphere  Mesoscale  Prediction 
System  (COAMPS).  COAMPS  features  a  nonhydrostatic  atmospheric  model  cou¬ 
pled  to  a  hydrostatic  ocean  model.  Initially,  COAMPS  is  intended  to  serve  as  a 
community  research  model  for  studying  various  mesoscale  phenomena  such  as  arc¬ 
tic  leads,  convection,  lake-eflfect  snowstorms  and  tropical  cyclones.  Some  of  these 
phenomena,  such  as  tropical  cyclones,  exchange  a  great  deal  of  energy  with  the 
ocean.  This  process  is  represented  within  COAMPS.  It  is  expected  that  the  study 
of  a  variety  of  mesoscale  phenomena  will  add  to  the  utility  of  COAMPS  as  it  trein- 
sitions  to  operational  use  over  the  next  few  years.  In  fact,  many  of  the  features  of 
the  COAMPS  forecast  models  will  eventually  transition  into  the  NORAPS  forecast 
model  in  our  coastal  data  assimilation  system. 

The  purpose  of  this  paper  is  to  describe  the  current  structure  of  COAMPS  and 
to  show  results  from  an  experiment  that  demonstrates  the  effect  of  ocean  changes 
on  the  structure  of  a  tropical  cyclone.  In  section  2  we  describe  the  nonhydrostatic 
atmospheric  model.  In  section  3  we  describe  the  hydrostatic  ocean  model.  The 
tropical  cyclone  experiments  with  COAMPS  are  described  in  section  4.  A  summary 
and  conclusions  are  given  in  section  5. 
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2.  Atmospheric  Model 


a.  Dynamics 

The  adiabatic  equations  are  developed  in  a  manner  similar  to  that  used  by 
Klemp  and  Wilhelmson  (1978).  Using  the  equation  of  state, 


P  —  pRdTv 

the  definition  of  the  virtual  temperature, 

r„  =  r(1.0  +  0.6089„) 


(1) 


(2) 


and  the  Exner  function, 


the  primitive  equations  can  be  rewritten  as: 


du 

a 


du  du  du  .  ^  _4 

=  -u- - v-^  -  xv-^  +  /u  +  I>„  +  Kb^*u 


dx 


(3) 

(4) 


dv 

dw 

'm 


:  dx’ 

..  dD. 

dv  dv 

V  Q  ■■  +  Rd  o  — 

dy  dy 

—  u; 

.  dx’ 

„  dDz 

dw 

dw 

dz  dz 

dv 


dz 


/u  +  +  Kb^*v  (5) 


dw 

’d^ 


rO’ 


+5(  J  +  0.608{9v  -  q„)  -qc-Qr-qi-  9^  (6) 


m  =  0 


(7) 


d9 

dt 


de  de 
""dx  '"dy 


dz 


where  Dz  is  the  three-dimensional  divergence, 

^  d(pO„u)  ,  d(p9vv)  ,  d(pd^.w) 


(8) 

(9) 
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p  is  the  pressure,  p  the  density,  the  g&s  constant  for  dry  air,  T  is  the  temperature, 
9v7  9c)  9r)  and  are  the  mixing  ratios  of  water  vapor,  cloud  droplets,  raindrops, 
ice  crystals  and  snowflakes,  respectively,  poo  is  a  constauit  reference  pressure,  Cp  is 
the  specific  heat  at  constant  pressure  for  the  atmosphere,  u,  v  and  w  are  the  wind 
components  in  the  x-,  y-  2uid  z-directions,  respectively,  /  is  the  Coriolis  parameter 
and  g  is  the  acceleration  due  to  gravity.  The  terms  £)„,  Dv,  D^,  and  De  represent 
subgrid  scale  mixing.  The  terms  involving  Kd  represent  divergence  damping  that 
control  the  amplification  of  sound  waves  (Skamarock  and  Klemp  1992)  where  Ku 
is  the  diffusion  coefficient.  Kb  is  a  coefficient  for  fourth-order  diffusion  used  to 
control  numerical  noise  when  the  equations  are  solved  on  a  discrete  grid.  The 
speed  of  sound  for  the  mean  state,  c,  is  defined  as: 


where  is  the  specific  heat  at  constant  volume  for  the  atmosphere.  In  eqs  (4)  - 
(8),  the  variables  k  and  6  have  been  decomposed  as: 

()  =  (■)  +  ()'  (11) 

where  the  overbar  denotes  the  initial  mean  state,  which  is  a  function  of  z  only, 
and  the  prime  denotes  deviations  from  this  mean.  The  mean  state  follows  the 
hydrostatic  relationship, 

—  -  __L- 

dZ  CpOy, 

The  equations  are  solved  on  a  staggered,  scheme  C  grid  ( Arakawa  and  Lamb,  1974). 
On  this  grid,  the  u-components  are  one-half  grid  interval  between  the  mass  points 
in  the  x-direction  and  the  u-components  axe  one-half  grid  interval  between  the 
mass  points  in  the  y-direction.  The  u7-component  is  staggered  in  a  similar  manner 
with  respect  to  the  mass  variables.  All  derivatives  are  computed  to  second-order 


5 


accuracy  except  the  horizontal  diffusion,  which  uses  approximations  of  fourth-order 
accuracy  to  damp  only  the  shortest  wavelengths.  High-frequency  time  oscillations 
are  controlled  using  the  Robert  (1966)  time  filter  using  a  smoothing  parameter  of 
0.2. 

b.  Subgrid  Scale  Mixing 

COAMPS  has  been  designed  to  use  horizontal  grid  spacing  ranging  from  an 
upper  limit  of  about  20  km  for  meso-;9  scale  applications  down  to  the  order  of  1  m 
for  Large  Eddy  Simulation  (LES)  applications.  At  the  higher  resolutions,  in  which 
the  horizontal  and  vertical  grid  spacings  are  typically  the  same,  the  subgrid  scale 
mixing  terms  must  represent  the  3-dimensional  turbulence  and  be  based  on  the 
volume- averaging  method  (Wyngaard  1982).  At  the  coarser  resolutions,  in  which 
the  horizontal  grid  spacing  can  be  an  order  of  magnitude  or  more  greater  than  the 
vertical  grid  spacing,  the  horizont2d  and  vertical  subgrid  scale  mixing  terms  must  be 
treated  separately  and  be  based  on  the  ensemble-averaging  method.  In  the  following 
sub-sections  we  describe  first,  the  volume- averaging  method  for  LES  applications; 
second,  the  ensemble- averaged  parameterization  for  meso-/?  applications  and;  third, 
the  background  diffusion  used  to  suppress  numerically  generated  noise. 

1)  LES  APPLICATIONS 

The  subgrid  scale  turbulence  terms.  Da,  for  any  variable  a,  are  represented 
eis; 

„  d(^)  ,  d(^)  .  d(^) 

*  ~dr 

The  fluxes  oV,  a'v'  and  a'w'  must  be  determined  to  close  the  system  of  equations. 
The  fluxes  are  written  in  terms  of  a  coefficient  of  eddy  viscosity,  K,  and  the  known 
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grid  point  variables.  The  momentum  fluxes  become, 


-7-7  „  /  dui  duj  \  2 

UiUj  -  -K^  +  —  j  +  -6. 

and  for  the  remaining  scalar  variables, 

da 


(14) 


u‘,a'  =  -Kt. 


dxi 


(15) 


where  Km  and  Kh  are  the  eddy  mixing  coefficients  for  momentum  and  heat,  re¬ 
spectively,  and  are  defined  as: 


Km  =  Smle^'^ 


(16) 


and 


Kk  = 


(17) 


e  is  the  subgrid  scale  turbulent  kinetic  energy  (tke)  following  Deardorff  (1980)  and 
a  represents  0,  q^,  qc^qT,qiiqs-  The  coefficient  Sm  is  set  to  a  constant  value  of  0.1, 
and 


and. 


(18) 


As  =  (AiAyAz)*/^ 


(19) 


where  Ai,  Ay  and  Az  are  the  grid  spacings  in  the  x-,  y-  and  z-directions.  The 
mixing  length,  /,  is  computed  as, 


/  =  0.76eM||^ 

9  oz 


-1/2 


(20) 


for  unsaturated  air  when  ||  >  0, 


'““■’''Kit 


-1/2 


(21) 
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for  saturated  air  when  ^  >  0  and 


1  =  As 


otherwise. 


The  turbulent  kinetic  energy,  e,  is  predicted  in  a  manner  similar  to  that 
described  by  DeardorfF  (1980).  Here, 

de  de  de  dt  ^  . 

where  BP  is  the  buoyancy  production  term,  SP  is  the  shear  production  term,  e 
is  the  dissipation  rate  and  all  the  other  variables  are  as  described  previously.  The 
buoyancy  production  is  defined  as: 

(24) 


for  nonsaturated  points  and 


BP  = 

9  dz 


for  saturated  points.  The  shear  production  is  defined  as: 


SP  =  (1^)  +(^)  +(^) 

\ox J  \dy )  \dz  / 

.  ( dv  / dw  du\^  (dw  duV 

\5x  ^  \dx  ^  dz )  ^  \  dy  ^  dz) 


-H 


du  dv  dw 
dx^  dy^  dz  ) 


and  the  dissipation  rate  is  defined  as: 


where 


C,  =  0.19  + 
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2)  MES0-;9  APPLICATIONS 

For  meso-/?  applications,  we  apply  eddy  viscosity  closure  separately  for  the 
horizontal  and  vertical  mixing  terms.  This  means  that  we  have  four  separate  eddy 
mixing  coefficients,  one  for  mixing  scalar  variables  in  the  horizontal,  KhHi  one  for 
mixing  momentum  in  the  horizontal,  Kmfj,  one  for  mixing  scalar  van'iables  in  the 
vertical,  Kf^v  and  one  for  mixing  momentum  in  the  vertical,  K^v-  The  eddy  mixing 
coefficients  in  the  vertical  are  as  defined  in  the  previous  section: 

Kkv  =  ShWt^^^  (29) 

and 


KmV  = 


(30) 


with  the  exception  that  /v,  represents  the  vertical  mixing  length,  defined  as: 

kz 


lv  = 


»0 


(31) 


where  k  is  the  von  Karman  constant  (0.4),  z  is  the  height  above  the  surf3w:e  and  Iq 
is  a  specified  constant.  The  prediction  of  e  is  performed  using  (23)  except  that  the 
shear  production  term  {SP)  is  simplified  as: 

2  /  n  \  ^  /  n  \2l 


SP  =  KmV 


(duV  (dvY  (dw\ 

+  (a;)  +  (^) 


(32) 


In  the  horizontal,  the  eddy  mixing  coefficients  are  based  on  the  deformation  of  the 
horizontal  velocity  field  (Smagorinsky  1963), 


KmH  — 


0.14/h^ 

V2 


fdu  ^  du 

^  dx)  dy) 


1/2 


(33) 


and 


KkH  —  ^KmH 


(34) 


9 


where  the  horizontal  mixing  length,  Ih  is  dehned  as; 


l„  =  (AxAy)'/*  (35) 

3)  BACKGROUND  DIFFUSION 

The  solution  of  continuous  equations  on  a  discrete  grid  inevitably  leads  to  the 
generation  of  spurious  gravity  waves  due  to  truncation  errors,  boundary  reflections, 
etc.  To  prevent  these  small  sources  of  energy  from  growing  and  contaminating  the 
forecasts,  we  have  included  a  background  fourth-order  diffusion  to  all  prognostic 
variables  except  for  tt'.  The  form  of  this  operator,  for  any  variable  a,  is  A'b  V^(V^q:), 
which  allows  us  to  compute  V^or  first,  use  this  on  the  grid  rows/columns  immedi¬ 
ately  adjacent  to  the  boundaries  and  then  compute  V^(V^a)  for  the  remaining  in¬ 
terior  points.  The  coeflRcient  is  chosen  such  that  the  ratio  (7<’BAta)/(AiAi^)=0.025 
for  the  second-order  points  and  (KsAto)/(Ax^At/^)=-0.0025  for  the  fourth-order 
points,  where  At,  is  the  time  step  of  the  model.  These  values  for  the  coefficient  have 
been  empirically  determined  so  that  they  produce  a  light  smoothing  of  the  shortest 
wavelengths  during  the  forecasts  with  minimal  effect  on  the  larger  wavelengths. 

c.  Moist  Physics 

Due  to  the  grid-spacing  used  in  COAMPS  (<20  km),  we  use  explicit  treat¬ 
ment  of  the  moist  physics,  rather  than  parameterizing  the  effects  of  convection.  The 
microphysical  processes  leading  to  conversion  between  water  vapor,  cloud  droplets, 
etc.  are  as  described  in  Rutledge  and  Hobbs  (1983)  and  therefore,  will  not  be  pre¬ 
sented  here.  The  inclusion  of  this  moist  physics  package  necessitates  equations  for 
water  vapor  (g„),  cloud  droplets  (qc),  raindrops  [q^],  ice  crystals  (q,)  and  snowflakes 
(qs)  in  addition  to  adding  a  latent  heating  term  (Qe)  to  the  thermodynamic  equa- 
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tion.  The  thermodynamic  and  moisture  equations  become: 


dt 


de  do  de  Qo  , 

=  - V- - w—  +  ~  +  De  +  Kb^  0 

ax  ay  oz  p 


(36) 


dqv 

dt 


—u 


dqv 

dx 


dqv  dqv  Sv  p.  I  i/'  VT"* 


(37) 


dqc 

dt 


— u 


dx 


—  V 


dqc 

dy 


+  D,^  +  Ks^^qc 


(38) 


dqr 

dt 


_  A  _  („  +  v,)^  +  4  +  +  KbV\, 

dy  '  dz  p  ^ 


dx 


(39) 


dqi 

dt 


dqi  dqi  dqi 

^  dx  '’dy  dz 


+  ~  +  D„  +  KB^Ui 

P 


(40) 


dqs 

dt 


-\-^  +  D„  +  KbV% 
P 


(41) 


where  the  terms  5«,  5e,  5r,  Si  and  5,  represent  sources  and  sinks  of  g„,  qc>  9r,  9« 


and  9,,  respectively,  Qg  represents  sources  and  sinks  of  heat  and  Vr  and  V,  are  the 


terminal  velocities  of  raindrops  and  snowflakes,  respectively. 


d.  Time-splitting 

The  complete  set  of  prognostic  equations  for  the  CO  AMPS  atmospheric  model 
are  (4)-(7),  (23)  and  (36)-(41).  Solving  these  fully-compressible,  nonhydrostatic 
equations  with  an  explicit  time  step  is  computationally  expensive  because  of  the 
speed  of  the  sound  waves.  However,  it  is  possible  to  isolate  the  terms  that  govern 
the  sound  waves  in  the  dynamiccil  equations  (i.e.,  the  left-hand-sides  of  (4)-(7)). 
Therefore,  a  short  time  step  can  be  used  for  the  left-hand-side  terms  while  a  longer 
time  step,  based  on  the  advective  speed,  can  be  used  for  the  right-hand-side  terms. 
Note  that  the  shorter  time  steps  are  not  necessary  for  the  thermodynamic,  moisture 
and  e  equations. 
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The  time-splitting  technique  used  here  follows  that  introduced  by  Klemp  and 
Wilhelmson  (1978).  For  a  small  time  step,  At*  and  a  large  time  step  At*,  (4)-(7) 
are  integrated  as: 


«— Ata+(n— l)ATa 


u 


-At„ 


(f-At«+{n-l)ATa) 


+Kd 
+RHS,.‘] 


(t) 


cl^V 

~dx  ) 

(t-Ata+(n-l)ATo) 


(42) 


^t-Ata+nAra 


V 


t— Ata+(n— l)ATa 


-At, 


I 


(t-Ata+(n-l)AT«) 


8  Oo  (*-^‘a+(n-l)AT«) 


(43) 


yjt-Ata+nAra 


W 


it— Ata+(n— r)ATo 


-At. 


(w) 


,\  (t-AtB+(n-l)ATo) 


+Kd 


(t-Ata+(n-l)ATa) 


dz 


(44) 


'y-Afo-^nATa  _  ^^/jt-Ato+(n-l)AT, 


^  /  r>  \l-Afa+nATo 

— 2  (^3) 


CppOv 


(45) 


where  RHS-u,  RHS„  and  RHS^,  are  the  right-hand-side  terms  in  (4),  (5)  and  (6),  re¬ 
spectively.  The  diffusion  coefficient  Ko  is  chosen  such  that  the  ratio  Ko^TalAxi  = 
0.015. 
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e.  Boundary  Conditions 


COAMPS  requires  the  use  of  lateral  as  well  as  upper  and  lower  boundary  con¬ 
ditions.  The  lateral  boundary  conditions  can  be  fixed,  periodic,  or  radiation  type 
similar  to  that  proposed  by  Orlanski  (1976).  When  using  the  radiation  bound¬ 
ary  conditions,  all  boundary  values  are  set  to  their  initial  values  at  inflow  points. 
At  outflow  boundary  points,  the  normal  velocity,  v„,  is  computed  using  upstream 
differencing  as: 


dv„  dvn 

IT  = 


(46) 


where  represents  the  velocity  normal  to  the  boundary  and  c.  is  some  estimate 
of  the  fastest  moving  gravity  wave  directed  out  of  the  domain.  In  COAMPS,  c, 
can  either  be  specified  or  computed  following  the  method  described  by  Miller  and 
Thorpe  (1981).  All  other  variables,  other  than  v„,  are  lineairly  extrapolated  to  the 
boundary  on  outflow.  It  should  be  noted  that  outflow  points  are  determined  by 
the  sign  of  (vn  +  c,).  The  upper  and  lower  boundaries  require  specification  of  the 
vertical  velocity  and  the  fluxes.  Currently,  we  set  tx;  =  0  at  both  the  upper  and 
lower  boundaries. 

The  lower  boundary  condition  on  the  vertical  mixing  terms  are: 


-(u'w')  =  «•  jyj 

(47) 

-(v'w')  =  “• 

(48) 

=  u,9. 

(49) 

-(^v)  =  w*9. 

(50) 

where  u,,  9,  and  qv»  are  velocity,  potential  temperature  and  mixing  ratio  scale 

values  and  |y|  =  (u^  -|-  at  the  first  grid  height  above  the  surface.  We  use  the 
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Louis  et  al.  (1982)  surface  flux  parameterization  to  compute  u.,  B,  and  q^,.  At 
the  upper  boundary,  the  normal  gradient  o^the  vertical  mixing  terms  is  assumed 
to  vanish  to  prevent  the  development  of  unrealistic  tendencies. 
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3.  Ocean  Model 


a.  Dynamics 

The  selection  of  the  ocean  model  for  COAMPS  was  determined  by  the  type 
of  experiments  we  anticipate  performing.  For  the  simulation  of  a  tropical  cyclone, 
Pudov  et  al.  (1978),  Hopkins  (1982),  Brooks  (1983)  and  Shay  and  Elsberry  (1987) 
all  demonstrated  that  the  oceanic  response  to  a  tropical  cyclone  can  extend  below 
the  thermocline,  hundreds  of  meters  below  the  surface.  Also,  proper  simulation 
of  the  response  of  the  ocean  mixed  layer  to  a  passing  tropical  cyclone  is  necessary. 
The  ocean  mixed  layer  exhibits  deepening  and  strong  mixing  while  the  SST  exhibits 
cooling  due  to  the  passage  of  a  tropical  cyclone  (Fisher,  1958;  Leiper,  1967;  Price, 
1981;  Shay  and  Elsberry,  1987;  Sanford  et  al.,  1987;  Brink,  1989).  Therefore, 
we  need  a  model  capable  of  representing  the  deep  (barotropic)  response  in  the 
ocean  as  well  as  resolving  the  mixing  processes  within  the  mixed  layer.  We  use  the 
incompressible,  hydrostatic  equations  similar  to  those  described  by  Chang  (1985) 
with  two  modifications.  First,  we  extend  the  axisymmetric  form  of  the  equations  to 
3-dimensions.  Second,  we  replace  the  predictive  equation  for  the  density  (p)  with 
predictive  equations  for  temperature  (T)  and  salinity  (5),  and  use  the  equation  of 
state  (Gill  1982)  to  diagnose  the  density  from  the  temperature,  salinity  and  pressure 
(p).  The  equations  become: 


dv 


dv  dv  r>A 

-u-T - V- - w-T - fu  +  Do.,  +  Kb^*v 

dx  dy  dz 


1  dpr 
Po  dy 


(53) 
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where  all  the  variables  are  defined  as  above  for  the  atmospheric  model  except  that 
here,  u,  v  and  w  are  ocean  current  speeds  in  the  x-,  y-  and  z-directions,  the  subscript 
5  refers  to  the  surface,  £?o„,  Uo„,  Dot  and  Dos  represent  subgrid  scale  mixing  of  the 
u-  and  v-components  and  the  temperature  and  salinity,  respectively,  and  are  defined 
below,  H  is  the  undisturbed  mean  height  of  the  ocean,  h  is  the  deviation  height 
from  this  mean,  referred  to  as  the  free  surface  height,  fp  symbolically  represents 
the  equation  of  state  and  pr  is  the  reference  pressure,  expressed  as: 

Pr{z)  =  p{z)  -  gpo{H  -  z)  (60) 

where  po  is  the  mean  ocean  density.  The  momentum  equations  have  been  decom¬ 
posed  into  a  vertical  mean  (barotropic  mode)  and  the  perturbation  from  this  mean 
(baroclinic  modes),  expressed  as, 

()  =  o+(r  (61) 
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Therefore,  the  mean  and  perturbation  reference  pressures  become: 


p,  =  —  +  gf/'z  +  p,gh  (62) 

and 

Pr  =  Pr  -  Pr  (63) 

respectively,  where  p,  is  the  surface  density. 

The  complete  set  of  equations  for  the  ocean  model  is  comprised  of  (51)-(58), 
(62)  and  (63).  As  in  the  atmospheric  model,  the  equations  are  solved  on  the 
staggered,  scheme  C  grid.  A  similar  staggering  is  used  in  the  vertical  with  tv 
defined  at  the  interface  between  the  layers  and  all  other  variables  are  defined  at  the 
middle  of  each  layer.  All  derivatives  are  computed  to  second-order  accuracy  except 
the  horizontal  diffusion,  which  uses  fourth-order  accurate  approximations  to  damp 
only  the  shortest  wavelengths.  High  frequency  time-oscillations  are  controlled  using 
the  Robert  (1966)  time  filter  using  a  smoothing  parameter  of  0.2. 

b.  Subgrid  Scale  Mixing 

The  subgrid  scale  mixing  in  the  COAMPS  ocean  model  is  similar  to  that 
used  in  the  FNOC  Thermodynamical  Oceem  Prediction  System  (TOPS)  (Clancy 
and  Martin,  1979)  and  is  based  on  the  level-2  turbulence  closure  of  Mellor  and 
Yamada  (1974).  The  mixing  vertical  terms  are  expressed  as: 

Dou  =  ^  (wV)  (64) 

Dov  =  ^  («^)  (65) 

= I  (^) 


(66) 


(67) 


D,s  =  3-  (=??) 


where, 


du  du 

%  " 

(68) 

dv 

%  ~  ~^‘^Tz 

(69) 

II 

1 

(70) 

dS  ..  dS 

Tz  = 

(71) 

where  /  is  the  mixing  length,  and  Sm  are  constants  computed  as  functions  of  the 
Richardson  number,  Koh  and  Kom  are  the  oce<in  eddy  mixing  coefficients  and 

q  =  (2«)‘'=  (72) 

where  e  is  the  turbulent  kinetic  energy.  The  details  for  the  calculations  of  5*  and  Sm 
Me  presented  in  Mellor  and  Yamada  (1974).  The  turbulent  kinetic  energy  is  com¬ 
puted  assuming  a  balance  between  the  shear  production,  the  buoyancy  production 
and  the  viscous  dissipation. 


”*  [[dzj  Vdz  j  j  V  dz )  151 


The  mixing  length,  I,  is  expressed  as: 


1  +  t 

where  k  is  the  von  Karman  constant  (0.4)  and  Iq  is  computed  as: 

,  0.1  Jl ft  \z\qdz 

f-H^dz 
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» 


where  H  is  the  depth  of  the  ocean. 

We  use  Iq  from  the  previous  time  step  when  solving  (73).  Optionally,  we  could 
iterate  each  time  step  in  the  calculations  of  q  and  /,  but  experiments  have  shown 
that  the  solutions  are  not  sensitive  to  this.  Also,  we  compute  the  mixing  terms  (68) 
-  (71)  as  adjustments  at  the  end  of  each  time  step  in  an  implicit  formulation  so 
that  we  are  not  constrained  by  the  stability  restrictions  on  Kgh  and  Kom  necessary 
when  using  an  explicit  formulation  of  these  terms.  The  horizontal  mixing  in  the 
ocean  model  consists  only  of  a  background  fourth-order  diffusion,  similar  to  that 
used  in  the  atmospheric  model  described  above. 


c.  Time-splitting 

One  of  the  complicating  factors  in  the  numerical  simulation  of  the  ocean  is 
the  wide  range  of  the  phase  speeds  of  the  various  modes  within  the  ocean.  Whereas 
the  external  (barotropic)  mode  is  on  the  order  or  300  m  s~^,  the  fastest  internal 
mode  is  on  the  order  of  10  m  s"^  and  the  advective  speeds  rarely  exceed  1-2  m  s“*. 
For  efficiency,  as  with  the  atmospheric  model,  we  use  a  time-splitting  technique 
to  reduce  the  number  of  computations.  We  solve  for  the  barotropic  mode  with 
a  small  time  step  (At^)  appropriate  for  the  barotropic  phcise  speed  and  a  larger 
time  step  (At^)  appropriate  for  the  phase  speed  of  the  fastest  internal  mode,  where 
A<o  is  an  integer  multiple  of  At,,.  Only  those  equations  that  govern  the  barotropic 
mode  are  integrated  with  the  smaller  time  step  (i.e.,  (51),  (53)  and  (57)).  Also,  the 
bracketed  terms  on  the  right-hand-side  of  (51)  and  (53)  only  need  to  be  evaluated 
every  large  time  step,  saving  additional  computer  time.  The  integration  of  u,  v  and 
h  are  carried  out  using  centered  differencing  in  time  and  space,  so  that, 


(^)‘-'^‘o-h(n+l)aTo 


(y)*-At<,+(n-l)ATo 


-2At„ 


Po 


t-Ato+(n-l)ATo1 


(76) 
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(y)t-Afo+(n+l)AT« 


(y)‘-^‘o+(n-l)AT„ 


-  2ATo 


Po 


t-Ato+(n-l)ATo1 


(77) 


ih) 


t-Ato+(n+l)ATo  _ 


(/l)*~A‘o+(n-l)ATo 

-2Ar 


dx 


t— Ato+nATo 


du  dv 


t— Ato+nATo 


(78) 


where  and  are  the  barred  and  bracketed  terms  on  the  right-hand-side  of  (51) 
and  (53),  respectively. 


d.  Boundary  Conditions 

The  lateral  boundary  conditions  in  the  ocean  model  are  similar  to  those  used 
in  the  atmospheric  model.  We  can  use  fixed,  periodic  or  radiation  type  boundary 
conditions.  When  using  the  radiative  boundary  conditions,  all  boundary  values  are 
set  to  their  initial  values  for  inflow  points.  At  outflow  points,  (46)  is  used  to  predict 
the  normal  velocity  component  while  all  other  variables  are  linearly  extrapolated 
to  the  boundary.  The  value  of  c.  is  fixed  at  10  m  s"^.  We  need  only  one  boundary 
condition  for  the  vertical  velocity,  since  we  allow  for  a  free  surface.  W’^e  set  u;=0  at 
the  lower  boundary  and  integration  of  the  continuity  equation  (59)  yields  w  at  all 
other  levels. 

The  upper  boundary  conditions  on  the  vertical  mixing  terms  are: 

-(W)  =  (79) 
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-(u;'5')  =  j  (u.q.5  -  P)  (82) 

where  T,  and  qy»  are  the  scale  values  of  velocity,  temperature  and  water  vapor 
computed  following  Louis  et  al.  (1982),  /)„  and  /)<,  are  the  densities  at  the  air- 
ocean  interface  of  the  atmosphere  and  ocean,  respectively.  Cop  is  the  specific  heat 
at  constant  pressure  for  the  ocean  and  P  is  the  precipitation  rate  at  the  surface.  At 
the  lower  boundary,  the  normal  gradient  of  the  vertical  mixing  terms  is  assumed  to 
vanish  to  prevent  the  development  of  unrealistic  tendencies. 
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4.  Model  Tests  and  Results 


a.  Background 

The  simulations  reported  here  are  based  on  data  from  Hurricane  Gilbert, 
which  moved  from  the  Caribbean  Sea  into  the  Gulf  of  Mexico  during  the  period  of 
8-19  September  1988.  During  its  existence,  Hurricane  Gilbert  established  the  record 
minimum  pressure  for  a  hurricane  at  888  mb  and  attained  a  700  mb  wind  speed  of 
89  m  s~^  (Willoughby  et  al.  1989).  Idealized  initial  conditions  are  used  to  con¬ 
struct  homogeneous  atmospheric  and  oceanic  structures,  typical  of  the  environment 
near  Gilbert,  to  demonstrate  whether  COAMPS  can  depict  the  development  of  the 
storm.  In  addition,  features  of  Gilbert,  such  as  concentric  eye  walls,  ais  discussed 
by  Black  and  Willoughby  (1992)  and  the  two-way  interaction  between  Gilbert  and 
the  ocean  are  examined. 

b.  Experiment  1:  Uncoupled  Simulation 

The  first  experiment  (EXPl)  performed  with  COAMPS  is  an  84  h  simulation 
of  Hurricane  Gilbert  using  a  fixed  sea  surfaw^e  temperature  (SST).  Therefore,  only 
the  atmospheric  model  grid  was  used  in  EXPl.  The  Cartesian  grid  consisted  of  61 
points  in  the  east-west  and  north-south  directions  with  a  horizontal  grid  spacing 
of  10  km  while  30  vertical  levels  were  used  with  a  constant  spacing  of  600  m.  The 
speed  of  sound  was  arbitrarily  set  to  100  m  s~^  in  the  model  so  that  the  large  time 
step  was  30  s  and  the  small  time  step  was  2.5  s.  Rauliative  boundary  conditions 
were  used  on  the  lateral  boundaries  with  c,  =  30  m  s~^.  Since  the  ratio  of  the 
horizontal  spacing  to  the  vertical  spacing  is  so  large  in  this  experiment,  the  meso- 
beta  form  of  the  subgrid  scale  mixing  was  employed.  Rather  that  using  (31)  to 
compute  the  vertical  mixing  length,  a  value  of  Iv  =  600  m  was  chosen  and  modified 
for  stable  conditions  based  on  (20)  and  (21).  The  vertical  mixing  length  was  further 
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constrained  so  that  it  was  no  less  than  12  m. 

The  environmental  data  was  chosen  to  try  to  simulate  the  rapid  deepening 
associated  with  Hurricaine  Gilbert.  A  climatological  SST  of  302.36  K  was  obtained 
from  the  Generalized  Digital  Environmental  Model  (GDEM,  Teague  et.  eil.,  1990) 
for  the  month  of  September  at  a  position  near  where  Gilbert  attained  it’s  maximum 
intensity  (lO.S^N,  83.9®W).  The  atmospheric  temperature  and  moisture  data  was 
obtained  from  the  Kingston,  Jamaica  sounding  valid  at  0000  UTC  September  11 
1988  (Fig.  1).  In  the  mid-  to  upper-levels,  the  relative  humidity  was  set  to  0.05 
per  cent  from  5500  m  to  7100  m  and  to  0.005  per  cent  above  7100  m.  A  Rankine 
vortex,  positioned  at  the  center  of  the  grid,  was  used  to  specify  the  initial  wind 
field.  Inside  the  radius  of  maximum  winds  (r,„ax)5  tangential  velocity  (vg)  was 
computed  as, 

=  W  (-^)  (l  -  (83) 

where  v^ax  is  the  maximum  tangential  velocity  (15  m  s“’),  r  is  the  distance  from 
the  center  of  the  storm  and  r^ox  is  set  to  90  km.  Outside  r^ax. 


where  x=0.6  is  an  empirical  constant  and  r,„  is  the  radius  of  influence  (240  km), 
beyond  which  vg  is  set  to  zero.  The  vertical  structure  of  the  tangential  wind  is 
given  by  (83)  and  (84)  below  10  km,  and  exponentially  decreases  to  zero  from  10 
km  to  13  km.  For  this  simulation,  no  mean  flow  w'as  used  so  that  the  entire  initial 
wind  field  was  comprised  solely  of  the  tangential  wind  from  the  Rankine  vortex. 

The  initial  pressure  field  was  computed  as  follows.  The  divergence  equation 
is  formed  from  the  horizontal  equations  of  motion  in  which  the  time- tendencies, 
diffusion  and  turbulent  mixing  terms  are  set  to  zero.  The  elliptic  equation  for  the 
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FIG.  1.  Skew  T-log  p  diagram  with  Kingston,  Jamaica  sounding  for  0000  UTC  11 
Sept  1988. 
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perturbation  pressure  becomes; 


Eq.  (85)  is  solved  for  each  model  level  assuming  that  =  0  on  the  lateral  bound¬ 
aries.  After  the  t'  fields  are  obtained,  they  are  adjusted  by  subtracting  the  value 
of  ir'  at  the  center  of  the  grid  for  each  level  from  all  points  on  that  level  to  insure 
that  the  data  nearest  the  center  of  the  grid  remains  as  close  to  the  original  speci¬ 
fied  values  as  possible.  The  potential  temperature  field  is  then  initialized  using  the 
vertical  equation  of  motion  (Eq.  6).  By  aissuming  that  uj,  the  time  tendency  of  w, 
diffusion  and  turbulent  mixing  are  zero,  (6)  can  be  written  ais: 

This  leads  to  a  maximum  adjustment  of  the  potential  temperature  of  about  -2°C 
between  10-13  km  above  the  surface  and  from  200  km  from  the  center  of  the  storm 
outward. 

The  initial  conditions  have  been  chosen  so  that  the  simulation  could  have 
the  best  chance  of  capturing  the  strong  intensification  associated  with  Gilbert. 
First,  no  environmental  shear  was  included  that  may  affect  the  vertical  structure 
of  the  storm.  Second,  the  SST  associated  with  the  position  of  maximum  intensity 
of  Gilbert  was  selected.  Third,  the  SST  was  not  allowed  to  vary  in  time  as  was 
observed  (Shay  et  al.,  1991)  which  could  lead  to  a  weakening  of  the  system.  This 
effect  will  he  studied  in  the  following  section.  Fourth,  the  model  grid  was  set  on 
an  f-plane  (17°N)  with  no  mean  wind  so  that  the  model  storm  stayed  symmetrical 
and  in  the  center  of  the  grid  to  minimize  effects  of  boundaries  and  asymmetries  in 
the  flow. 

The  84  h  COAMPS  forecast  minimum  pressure  and  the  maximum  wind  speeds 
for  EXPl  and  the  observed  minimum  pressure  for  Gilbert  are  shown  in  Fig.  2. 
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FIG.  2.  a)  Observed  central  pressure  (mb)  from  0000  UTC  11  Sept  1988  to  1200 
UTC  14  Sept  1988  for  Hurricane  Gilbert  (solid)  and  COAMPS  central  pressure 
forecaist  for  EXPl  (dash).  Values  are  plotted  in  6  h  increments,  (b)  As  in  (a) 
except  for  maximum  wind  speed  (m  s~^). 


The  initial  time  of  the  model  simulation  is  taken  to  be  0000  UTC  11  September 
1988  since  this  is  the  time  of  the  sounding  used  for  the  initial  conditions.  At  this 
time,  Gilbert  was  observed  to  have  a  minimum  pressure  of  approximately  983  mb. 
Over  the  next  36  h,  Gilbert  slowly  deepened  to  about  960  mb,  then  remained  at 
approximately  this  intensity  for  the  next  12  h  as  it  moved  over  the  island  of  Jamaica. 
Gilbert  deepened  rapidly  over  the  next  24  h  to  attain  a  record  low  pressure  of  888 
mb  at  0000  UTC  14  September,  corresponding  to  72  h  in  the  model  simulation.  At 
this  time,  Gilbert  attained  a  maximum  wind  speed  of  approximately  83  m  s~^. 

Two  distinct  stages  of  development  of  Gilbert  in  EXPl  can  be  seen  from 
Fig.  2.  In  the  first  stage,  defined  as  the  first  24  h  of  the  simulation,  the  central 
pressure  remained  nearly  constant  at  1000  mb.  During  this  stage,  the  model  storm 
developed  low-level  inflow  in  response  to  the  surface  friction  and  a  warm,  moist 
PEL  due  to  the  sensible  and  latent  heat  fluxes.  In  the  second  stage,  defined  as  the 
last  60  h  period  from  24  h  to  84  h  of  the  simulation,  the  model  storm  exhibited 
rapid  deepening  from  1000  mb  to  the  minimum  of  874  mb  at  81.5  hours.  Although 
the  pressure  fell  during  this  second  stage,  there  were  periods  of  pressure  rises  and 
quasi-steady  state  behavior  imbedded  within  the  general  pattern  of  pressure  falls. 
A  description  of  these  two  stages  of  the  model  storm  development  follows. 

During  the  first  stage,  the  model  storm  developed  inflow  in  the  lowest  levels 
through  surface  friction.  This  led  to  an  initial  spin-down  of  the  storm  circulation  for 
approximately  the  first  12  h  until  the  induced  upward  motion  initiated  convection. 
This  convection  formed  near  the  radius  of  maximum  winds,  initially  set  at  90  km, 
but  rapidly  contracted  inward.  East-west  cross-sections  of  the  u— component  (i.e., 
tangential  wind)  for  0  h,  12  h  and  24  h  in  Fig.  3  show  the  initial  decay  of  the 
circulation  followed  by  a  spin-up  with  the  maximum  tangential  velocity  over  30  m 
s~^  at  24  h  only  40  km  from  the  center  of  the  circulation. 
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FIG.  3.  East-west  vertical  cross-sections  of  the  tangential  velocity  (m  s~^)  taken 
through  the  center  of  the  grid  for  EXPl  for  a)  0  h,  b)  12h,  and  c)  24h.  Contour 
interval  is  5  m  s~^ 
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FIG.  3,  continued. 
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The  maximum  wind  speed  in  the  model  started  increasing  at  19  h  of  the  sim¬ 
ulation  during  which  time  the  minimum  pressure  remained  constant.  The  increase 
of  the  maximum  wind  prior  to  the  central  pressure  fall  is  similar  to  that  reported 
by  Rosenthal  (1977).  The  increase  in  wind  speed  occurs  about  12  h  after  the  first 
convective  towers  formed  in  the  eye  wall.  However,  it  is  only  2  hours  after  the 
development  of  ice  crystals  atop  the  eye  wall  convection.  The  added  heat  release 
due  to  the  formation  of  the  ice  clouds  appears  to  be  a  significant  factor  in  causing 
the  deepening  of  the  system  as  also  found  by  Tripoli  (1992b). 

The  second  stage  of  development  of  EXPl  was  characterized  by  the  rapid 
fall  of  the  central  pressure  from  approximately  1000  mb  to  874  mb.  During  this 
period,  the  model  storm  developed  many  features  that  have  been  observed  in  trop¬ 
ical  cyclones.  Model  generated  radar  reflectivities  at  3300  m,  computed  following 
Rutledge  and  Hobbs  (1983)  formulation,  for  every  12  h  from  24  h  to  84  h  (Fig.  4) 
illustrate  the  structure  of  the  convection  in  the  model  storm.  Strong  convective 
cells  (greater  than  40  dBz)  are  evident  in  a  ring  surrounding  the  eye  at  a  radius  of 
20-40  km.  No  precipitation  is  found  in  the  eye.  By  48  h,  the  appearance  of  spiral 
bands  of  convection  at  radii  ranging  from  50-200  km  is  noted.  Within  these  bands, 
convective  cells  propagate  inward  with  reflectivities  often  equal  to  cells  within  the 
eye  wall.  The  pattern  of  4  symmetric  spiral  bands  is  due  to  the  use  of  an  f-plane, 
no  mean  flow  and  the  storm  being  located  in  the  center  of  a  square  grid. 

The  cross-section  of  the  vertical  velocity  at  84  h  (Fig.  5)  illustrates  the  upward 
motion  associated  with  the  eye  wall  convection.  In  general,  the  maximum  upward 
motion  within  the  eye  wall  varied  in  time  from  2-6  m  s“*.  This  was  also  found  to  be 
the  case  for  the  individual  convective  elements  within  the  spiral  bands.  Downward 
motion  within  the  eye  exhibited  a  maximum  of  about  -0.5  m  s“^  near  9  km  with 
weaker  subsidence  extending  nearly  to  the  surface.  The  vertical  motion  associated 
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FIG.  4,  continued. 
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FIG.  4,  continued. 
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FIG.  4,  continued. 
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FIG.  4,  continued. 
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FIG.  5.  East- west  vertical  cross-sections  of  the  vertical  velocity  (m  s“^)  taken 
through  the  center  of  the  grid  for  EXPl  at  84  h.  Contour  interval  is  1  m  s~^. 
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with  the  eye  wall  convection  slanted  outward  with  height  as  has  been  found  by 
Rosenth2d  (1977). 

The  radial  wind  component  at  48  h  is  shown  in  Fig.  6a.  Inflow  of  over  10 
m  s~^  dominates  below  1  km  near  the  eye  wall  indicating  the  large  convergence  in 
this  region.  A  secondary  maximum  in  low-level  convergence  is  found  at  about  100 
km  from  the  center  of  the  storm  in  conjunction  with  the  spiral  bands.  In  general, 
5  m  s~^  inflow  exists  throughout  the  lowest  kilometer.  Outflow  is  found  at  10-13 
km  with  speeds  in  excess  of  10  m  s~'.  By  84  h,  a  much  stronger  radial  circulation 
has  developed  (Fig.  6b).  The  near-surface  inflow  is  now  over  20  m  s”^  with  inflow 
over  5  m  s“^  up  to  nearly  4  km.  A  strong  outflow  of  over  25  m  s“^  exists  at  about 
12  km. 

It  was  shown  earlier  that  the  central  pressure  of  the  model  storm  increased 
somewhat  after  about  48  h  of  simulation.  This  coincides  with  the  formation  of 
the  spiral  bands  of  convection.  These  bands  acted  as  a  dynamic  part  of  the  sys¬ 
tem,  occasionally  contracting  inward  toward  the  eye  with  noticeable  changes  on 
the  central  pressure  and  maximum  wind  speed.  Willoughby  et  al.  (1982)  noted 
that  concentric  eye  walls  have  been  observed  in  intense,  symmetric  hurricanes.  Fur¬ 
thermore,  the  outermost  concentric  ring  of  convection  typically  contracts  inward, 
thereby  weakening  the  innermost  ring  of  convection  until  the  outer  ring  replaces  the 
inner  ring  as  the  eye  wall.  During  the  period  of  the  contraction  of  the  outer  ring,  the 
storms  were  observed  to  weaken  in  intensity,  with  rapid  intensification  following  the 
replacement  of  the  eye  wall  by  the  outer  ring.  Other  observational  studies  (Fortner 
1958,  Jordan  and  Schatzle  1961,  Jordan  1966,  Holliday  1977  and  Willoughby  1992) 
show  that  the  eye  of  the  storm  is  somewhat  larger  after  the  replacement  of  the  inner 
eye  wall  by  the  outer  eye  wall. 

Shapiro  and  Willoughby  (1982)  demonstrated  that  a  heat  source,  such  from 
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FIG.  6.  East- west  vertical  cross- sections  of  the  radial  wind  (m  s"^)  taken  through 
the  center  of  the  grid  for  EXPl  at  (a)  48  h,  and  (b)  84  h.  Contour  interval  is  5  m 
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convection,  near  a  maximum  in  tangential  wind  leads  to  larger  heat  falls  inside  the 
radius  of  maximum  winds  than  outside.  This  leads  to  a  propagation  inward  of  the 
tangential  wind  maximum  and  the  associated  heat  source.  The  contraction  of  the 
wind  maximum  often  leads  to  a  destruction  of  an  inner  wind  maximum  and  an 
eventual  replacement  by  the  outer  wind  maximum. 

Several  occurrences  of  concentric  eye  walls  were  present  in  EXPl  after  48  h 
of  the  simulation.  One  of  the  most  significant  occurred  during  the  period  from  48 
h  to  60  h.  Model  generated  radar  reflectivities  for  the  3300  m  height  at  48  h,  50  h, 
52  h,  54  h,  56  h  and  57  h  are  given  in  Fig.  7.  At  48  h,  a  ring  of  convection  is  found 
around  the  eye  with  the  newly  formed  spiral  bands  over  100  km  from  the  center 
of  the  storm.  The  maximum  reflectivity  in  the  eye  wall  is  found  20  -  30  km  from 
the  eye.  Over  the  next  6  hours,  the  spiral  bands  slowly  converge  in  toward  the  eye, 
forming  a  nearly  closed  ring  of  convection  around  the  eye  at  a  radius  of  about  70 
km  by  54  h.  By  this  time,  the  radar  reflectivities  are  larger  in  the  outer  eye  waU 
than  in  the  inner  eye  wall.  At  56  h,  it  appears  that  the  outer  and  inner  eye  walls 
are  merging  into  one.  Individual  cloud  elements  make  up  a  relatively  disorderly 
pattern  of  convection  from  20-60  km  from  the  eye.  One  hour  later,  at  57  h,  the 
concentric  rings  have  merged  into  one  eye  wall,  with  the  maximum  reflectivities 
now  about  35-40  km  from  the  eye. 

The  central  pressure  and  maximum  wind  speeds  for  the  time  period  of  48 
h  to  60  h  are  shown  in  Fig.  8  with  a  higher  temporal  resolution  than  weis  used 
in  Fig.  2a.  From  48  h  to  50.5  h,  the  period  in  which  the  spiral  band  was  still 
at  least  80  km  from  the  eye,  the  pressure  continued  to  fall  and  the  wind  speed 
continued  to  rise.  From  50.5  h  to  54.25  h,  the  time  period  over  which  the  outer 
eye  wall  contracted  from  80  km  to  60  km,,  the  central  pressure  rose  6.5  mb  and 
the  maximum  wind  decreased  8.5  m  s~^.  Only  after  54.25  h,  when  the  outer  and 
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FIG.  7.  Model  generated  radar  reflectivites  at  3300  m  for  EX  PI  over  central  21  x 
21  grid  points  at  (a)  48  h,  (b)  50  h,  (c)  52  h,  (d)  54  h,  (e)  56  h,  and  (f)  57  h. 
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FIG.  7,  continued. 
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FIG.  7,  continued. 
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FIG.  7,  continued. 
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FIG.  8.  (a)  Central  pressure  (mb)  forecast  for  EXPl  from  48-60  h.  Values  plotted 
every  time  step  (30  s).  (b)  As  in  (a)  except  for  maximum  wind  speed. 
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inner  eye  walls  began  to  merge,  did  the  central  pressure  begin  to  fall  agadn  with 
a  corresponding  increase  in  the  maximum  winds.  Rapid  intensification  occurred 
during  the  one  hour  time  period  from  57  h  to  58  h  with  the  central  pressure  falling 
9.6  mb  from  924.3  mb  to  914.7  mb  and  the  maximum  winds  increasing  from  69.3 
m  s"^  to  74.0  m  s~^.  This  intensification  was  short-lived,  however,  when  the  spiral 
bands  st2U'ted  contracting  inward  toward  the  eye  again  by  60  h.  As  noted  above, 
the  radius  of  the  eye  increased  from  about  25  km  to  40  km  during  the  merging  of 
the  two  concentric  eye  walls. 

An  examination  of  the  overall  deepening  of  the  model  storm  highlights  the 
overall  effects  of  spiral  bands  on  the  storm  development.  It  was  noted  above  that 
the  contraction  of  the  spiral  bands  into  an  outer  concentric  ring  that  contracted 
inward  resulted  in  a  temporary  weakening  of  the  system.  A  re-examination  of  Fig. 
6  shows  that  overall  inflow  and  outflow  of  the  storm  is  controlled  significantly  by 
the  presence  of  the  outer  spiral  bands  at  48  h.  However,  at  84  h,  when  the  system 
has  just  undergone  a  significant  deepening  of  over  20  mb  during  the  previous  12  h, 
there  is  noticeable  decrease  in  spiral  band  convection  (Fig.  7).  The  cross-sections  of 
the  radial  components  at  48  h  and  84  h  (Fig.  6)  illustrate  the  effect  that  the  spiral 
bands  have  on  the  storm  circulation.  At  48  h,  the  inflow  is  impeded  into  the  eye  by 
the  convergence  near  100  km  radius.  In  addition,  the  strong  outflow  aloft  appears  to 
be  associated  most  closely  with  the  outer  band  of  convection  rather  than  the  inner 
eye  wall.  It  appears  that  the  proximity  of  the  outer  band  of  convection  is  acting 
to  retard  the  development  and/or  maintenance  of  the  inner  eye  wall.  However,  at 
84  h,  without  the  presence  of  a  significant  outer  eye  wall,  the  low-level  inflow  is 
considerably  stronger  and  the  upper-level  outflow  is  tied  to  the  eye  wall  convection 
and  proceeds  undisturbed  out  to  the  edge  of  the  forecast  grid. 
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c.  Experiment  2:  Coupled  Simulation 

The  second  experiment  (EXP2)  performed  with  COAMPS  is  an  84  h  coupled 
simulation  of  hurricane  Gilbert  using  both  the  atmospheric  and  oceanic  models. 
The  atmospheric  model  was  set  up  exactly  as  in  EXPl.  The  ocean  model  also 
used  a  61  X  61  Cartesian  grid  so  that  the  grid  points  for  the  atmosphere  and  ocean 
were  coincident.  The  ocean  model  used  20  verticad  levels,  with  the  spacing  set  at 
5  m  for  the  uppermost  layer  and  each  successively  lower  layer  being  1  m  more  in 
depth.  Therefore,  the  total  ocean  depth  was  290  m.  The  large  time  step  for  the 
ocean  model  was  set  to  1200  s  and  the  small  time  step  was  30  s.  The  coupling  was 
achieved  by  integrating  the  atmospheric  model  first  for  40  time  steps,  i.e.,  for  one 
ocean  time  step.  The  average  surface  fluxes  of  heat,  momentum  and  moisture,  and 
the  precipitation  over  these  40  time  steps  were  then  used  as  the  upper  boundary 
conditions  for  the  ocean  subgrid  scale  mixing  processes. 

The  atmosphere  was  initialized  with  a  Rankine  vortex  exactly  as  described 
above  for  EXPl.  The  ocean  model  was  assumed  to  be  at  rest  initially,  and  the 
GDEM  profiles  of  temperature  and  salinity  (Fig.  9)  that  were  used  to  obtain  the 
SST  in  EXPl  were  used  in  EXP2  to  define  the  homogeneous  vertical  thermody¬ 
namic  structure  of  the  ocean.  Therefore,  the  starting  values  for  the  atmospheric 
forecast  in  EXP2  were  exactly  those  used  in  EXPl.  However,  in  EXP2,  the  SST 
(i.e.,  the  temperature  of  the  top  layer  of  the  ocean)  wjis  allowed  to  vary,  with  the 
changes  being  due  to  the  influence  of  the  surface  energy  exchanges  and  precipita¬ 
tion. 

The  COAMPS  forecasts  for  the  central  pressure  and  the  maximum  wind  speed 
for  EXPl  and  EXP2  are  presented  in  Fig.  10.  As  in  EXPl,  the  coupled  simulation 
requires  about  24  h  for  the  initial  spin-up  before  the  hurricane  starts  to  deepen. 
Through  the  first  60  h  of  the  simulations,  EXP2  tends  to  be  somewhat  weaker, 
with  a  central  pressure  difference  between  the  two  simulations  of  generally  5-10  mb 
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FIG.  9.  (a)  Initial  ocean  temperature  (K)  profile  for  EXP2.  (b)  As  in  (a)  except 
for  salinity  (ppt). 
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FIG.  10.  (a)  COAMPS  central  pressure  (mb)  forecast  for  EXPl  (solid)  and  EXP2 
(dash),  (b)  As  in  (a)  except  for  maximum  wind  speed  (m  s“^). 
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and  a  maximum  wind  speed  difference  of  2-10  m  from  24  h  to  60  h.  After  60 
h,  more  significant  differences  are  found.  While  in  the  fixed  SST  experiment  the 
central  pressure  dropped  from  912  mb  to  874  mb  after  60  h,  in  EXP2  the  central 
pressure  rose  from  920  mb  to  936  mb  at  72  h  before  it  intensified  somewhat  to 
930  mb  by  84  h.  At  84  h,  the  difference  in  central  pressure  between  EXPl  and 
EXP2  is  53.5  mb.  Similar  tendencies  are  found  in  the  maximum  wind  speeds  with 
a  difference  of  21  m  s“^  at  84  h. 

The  reason  for  the  behavior  of  EXP2  relative  to  EXPl  can  be  seen  by  looking 
at  the  SST  change  near  the  center  of  the  storm  (at  10  km  radius)  in  EXP2  (Fig. 
11 ).  During  the  first  42  h,  the  SST  exhibited  cooling  of  less  than  0.5°C.  After  42  h, 
the  cooling  increased  dramatically,  with  SST  changes  of  -5.7°C  by  60  h  and  -8.0“C 
at  84  h.  At  outer  radii,  the  SST  cooling  was  generally  no  less  than  -0.5“C.  From 
these  SST  changes,  it  is  understandable  that  there  are  relatively  small  differences 
between  the  strength  of  the  hurricane  in  EXPl  and  EXP2  during  the  first  42  h. 
However,  since  the  SST  dropped  over  5°C  from  42  h  to  60  h  with  no  significant 
weakening  occurring  in  the  strength  of  the  hurricane  in  EXP2  until  after  60  h,  it 
appears  that  there  is  a  time  lag  on  the  order  of  12  h  between  strong  surface  cooling 
near  the  eye  wall  and  corresponding  changes  in  the  storm  intensity.  This  result  is 
similar  to  that  found  by  Anthes  and  Chang  (1978). 

The  distribution  of  the  SST  changes  for  EXP2  over  the  entire  domain  at  84  h 
is  given  in  Fig.  12.  The  maximum  change  of  8.3°C  is  found  about  20  km  from  the 
storm  center.  Cooling  of  up  to  3‘’C  extends  to  a  radius  of  about  80  km  while  much 
of  the  remainder  of  the  domain  exhibits  an  SST  reduction  of  generally  about  0.5°C. 
Some  enhanced  cooling  (of  up  to  1.2°C)  is  noted  under  the  location  of  the  spiral 
bands.  This  can  be  traced  to  the  cooling  of  the  near-surface  air  temperature  from 
the  evaporation  of  raindrops,  resulting  in  a  slight  increase  in  the  surface  sensible 
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FIG.  11.  SST  change  (K)  from  0-84  h  10  km  from  center  of  the  storm  in  EXP2. 
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FIG.  12.  SST  change  (K)  at  84  h  for  the  full  forecast  domain  in  EXF2. 


heat  flux  (not  shown). 

Cross-sections  illustrating  the  time  evolution  of  the  ocean  temperature  change 
are  given  in  Fig.  13.  At  48  h,  the  maximum  cooling  of  about  4“C  is  found  at  a 
depth  of  62  m,  near  the  initial  position  of  the  thermocline.  The  cooling  proceeds 
so  that  at  84  h,  the  maximum  cooling  of  -9.5°C  is  also  found  at  a  depth  of  62  m 
with  nearly  all  changes  of  more  than  2.0“  C  confined  to  within  50  km  of  the  center 
of  the  storm.  At  the  surface,  the  maximum  cooling  of  -8.3“C  is  found  20  km  from 
the  storm  center  with  -7.3“ C  cooling  directly  under  the  eye  of  the  hurricane. 

The  ocean  current  speeds  and  direction  induced  by  the  hurricane  at  84  h  are 
shown  in  Fig.  14.  Maximum  speeds  of  3.1  m  are  found  30  km  from  the  eye, 
directly  under  the  surface  maximum  winds.  A  strong  pattern  of  divergence  exists 
within  the  cyclonic  circulation  within  a  50  km  radius  of  the  eye.  This  divergence 
leads  to  upwelling  in  the  ocean,  with  maximum  vertical  velocities  as  high  as  0.32 
cm  s”^  This  is  about  twice  as  large  as  the  value  found  by  Chang  and  Anthes  (1979) 
for  a  storm  moving  at  5  m  s“^.  Rather  than  occurring  directly  under  the  eye,  the 
strongest  vertical  velocity  occurred  at  a  radius  of  10-20  km  with  weaker  upward 
motion  under  the  eye  (not  shown). 

The  salinity  of  the  ocean  at  84  h  in  the  uppermost  layer  of  the  ocean  is  given 
in  Fig.  15.  Inside  a  radius  of  about  50  km,  the  salinity  has  increased  from  an  initial 
value  of  36.0  ppt  to  36.6  ppt.  This  is  a  result  of  the  combination  of  vertical  mixing 
and  upwelling  bringing  larger  sub-surface  values  to  the  surface  and  large  latent  heat 
fluxes  (>  3000  W  m”^)  inside  the  radius  of  maximum  winds  resulting  in  a  loss  of 
fresh  water  from  the  top  ocean  model  layer.  Another  area  of  change  is  found  under 
the  spiral  bands  where  the  salinity  decreeised  from  36.0  ppt  to  35.1  ppt  over  the  84 
h  of  the  simulation.  It  was  shown  earlier  that  these  areas  also  exhibited  a  decrease 
in  SST.  The  decrease  in  salinity  here  is  due  to  the  imbalance  between  the  relatively 
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FIG.  13.  East-west  vertical  cross-sections  of  the  ocean  temperature  change  (K)  for 
EXP2  at  (a)  48  h  and  (b)  84  h.  Contour  interval  is  1  K. 
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FIG.  15.  84  h  EXP2  forcca.st  of  salinity  (ppt)  at  2.5  m  depth. 


weak  latent  heat  flux  in  this  area  {  250  W  m“*)  and  the  total  precipitation  over  84 
h. 

The  SST  changes  of  over  8®C  reported  here  are  somewhat  larger  than  the 
5-6®C  reduction  actually  observed  for  Hurricane  Gilbert  (Shay  et  al.  1991).  This 
can  be  partially  explained  by  the  fact  that  in  EXP2,  the  model  storm  remained 
stationary  for  all  84  h  of  the  simulation  and  therefore,  adl  the  mixing  and  upwelling 
responsible  for  the  SST  reduction  remaned  concentrated  over  the  same  area.  In 
reality,  Gilbert  moved  westward  to  west-northwestward  at  a  speed  of  about  7  m 
s“^  and  the  SST  cooling  processes  were  isolated  for  a  given  area  for  a  much  shorter 
period  of  time.  These  findings  indicate  that  SST  cooling  may  have  little  effect  on 
hurricane  strength  for  storms  with  a  movement  of  as  little  as  1-2  m  s~*.  Since  the 
storm  is  moving  it  will  continuously  be  entering  areas  of  undisturbed  SST,  unless 
of  course,  the  storm  has  looped  back  over  a  previous  position  or  another  storm’s 
track.  Also,  since  there  appears  to  be  a  time  lag  of  perhaps  12  h  between  the  strong 
SST  cooling  and  a  response  to  the  hurricane  strength,  a  forward  speed  of  1-2  m  s~* 
will  enable  the  eye  wall  of  the  storm  to  move  away  from  the  cooled  area  before  it 
can  have  an  impact  on  the  eye  wall. 

An  examination  of  model-generated  radar  reflectivities  for  EXP2  (Fig.  16) 
compared  to  those  described  earlier  for  EXPl  (Fig.  4)  illustrate  the  differences 
in  the  structure  of  the  hurricane  between  the  two  simulations.  In  both  EXPl 
and  EXP2,  similar  patterns  are  observed  up  to  48  h.  The  convection,  initially 
at  radii  from  50-100  km,  collapses  inward  to  form  an  eye  wall  with  maximum 
reflectivities  within  20-40  km  of  the  center  of  the  storm.  Spiral  bands  are  evident 
in  both  simulations  by  48  h,  somewhat  stronger  in  EXPl,  presumably  since  the 
SST  has  decreased  as  much  as  0.5°C  in  EXP2  over  much  of  the  domain.  After  48 
h,  significant  differences  are  noted  with  the  eye  wall  structure.  First,  at  60  h  in 
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FIG.  16.  As  in  Fig.  4  except  for  EXP2. 
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FIG.  16,  continued. 
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FIG.  16,  continued. 


EXP2,  it  is  evident  that  the  eye  wall  is  starting  to  expand  in  response  to  the  SST 
cooling.  Also,  the  strength  of  the  echo  at  radii  from  10-20  km  is  weaker  in  EXP2. 
At  72  h  and  84  h,  the  echo-free  region  associated  with  the  eye  is  noticeably  larger 
in  EXP2.  At  84  h,  at  the  time  of  the  maximum  SST  reduction,  EXP2  exhibits 
eye  wall  that  is  thinner  and  weaker  than  that  in  EXPl.  In  spite  of  the  strong  SST 
cooling  near  the  center  of  the  circulation,  the  overall  structure  of  EXP2,  outside  a 
radius  of  approximately  50  km,  remains  the  same  with  the  spiral  bands  continuing 
to  feed  convective  cells  inward  toward  the  eye  wall.  The  significant  change  due  to 
the  SST  cooling  is  with  the  strength  and  position  of  the  eye  wall. 

The  differences  in  the  EXPl  and  EXP2  tangential  velocities  at  48  and  84 
h  can  be  seen  by  examining  Fig.  17.  At  48  h,  both  EXPl  and  EXP2  exhibit 
maximum  tangential  velocities  of  62  m  s“^  at  a  radius  of  20  km.  The  tangential 
wind  profiles  are  nearly  constant  with  height  up  to  5  km  with  the  winds  decreasing 
above  this  level.  By  84  h,  the  most  significant  differences  are  found  in  the  eye 
wall  region.  Besides  EXP2  exhibiting  a  maximum  tangential  wind  of  67  m  s~^  as 
opposed  to  85  m  s~^  in  EXPl,  the  radius  of  the  maximum  wind  in  EXP2  is  10  km 
larger  than  in  EXPl.  Also,  the  gradient  of  the  tangential  wind  from  the  maximum 
to  the  eye  is  significantly  weaker  in  EXP2.  At  radii  beyond  about  50  km,  only 
small  differences  exist  in  the  tangential  wind  profiles.  This  indicates  that  while  the 
intensities  (based  on  the  maximum  wind  speed)  of  the  EXPl  and  EXP2  storms  are 
different,  their  strengths  (based  on  winds  outside  the  radius  of  maximum  winds) 
appear  to  be  quite  similar.  This  suggests  that  modification  of  the  SST  by  a  tropical 
cyclone  may  have  little  or  no  effect  on  the  subsequent  track  of  the  storm  due  to 
any  forced  changes  in  structure  since  the  ‘^ST  changes  will  only  affect  the  intensity 
and  not  the  overall  strength  of  the  storm  (Holland  1984,  DeMaria  1985). 
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FIG.  17.  East-west  vertical  cross-sections  of  tangential  velocity  (m  s~^)  taken  from 
the  left  boundary  to  the  center  of  the  grid  for  (a)  EXPl  at  48  h,  (b)  EXP2  at  48 
h,  (c)  EXPl  at  84  h,  and  (d)  EXP2  at  84  h.  Contour  interval  is  10  m  s~^. 
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FIG.  17,  continued. 
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5.  Conclusion 


A  nonhydrostatic,  compressible  atmospheric  model  has  been  coupled  to  a 
hydrostatic  ocean  model.  This  coupled  ocean/atmosphere  mesoscale  prediction 
system  (COAMPS)  is  flexible  in  that  it  is  written  in  standard  FORTRAN  77  code 
and  the  number  of  grid  points  and  the  resolution  of  the  grid  can  be  set  by  the 
user  for  all  three  dimensions.  The  models  in  COAMPS  can  be  run  simultaneously 
for  fully  coupled  simulations  or  either  one  can  be  run  on  a  stand-alone  basis  for 
an  atmospheric  or  oceanic  simulation.  Coupling  of  the  fluxes  across  the  air-water 
interface  is  accomplished  by  averaging  the  atmospheric  fluxes  for  all  time  steps 
within  one  ocean  time  step.  The  coupling  includes  exchanges  of  heat,  momentum 
and  moisture  plus  the  effects  of  precipitation  on  ocean  salinity.  COAMPS  currently 
does  not  include  radiation,  the  effects  of  the  precipitation  temperature  on  the  ocean 
temperature,  the  evaporation  of  ocean  spray,  or  the  effect  of  ocean  waves  on  the 
surface  drag. 

The  atmospheric  model  has  been  designed  for  mesoscale  and  large-scale  eddy 
simulations.  An  explicit  moist  physics  routine  allows  for  the  prediction  of  water 
vapor,  cloud  water,  raindrops,  ice  crystals  and  snowflakes.  A  1-1/2  order  subgrid 
scale  closure  model  is  used  to  predict  the  turbulent  kinetic  energy.  A  time-splitting 
technique  is  used  for  efficiency  in  handling  the  sound  and  gravity  waves. 

The  ocean  model  uses  a  time-splitting  scheme  to  separate  the  barotropic  and 
baroclinic  terms.  Predictive  equations  for  u-,  v-,  temperature,  salinity  and  the 
height  of  the  free  surface  are  included.  The  model  uses  a  Richardson  number  based 
diagnostic  equation  for  the  subgrid  scale  mixing  terms. 

Two  idealized  case  studies  of  Hurricane  Gilbert  have  been  performed  using 
COAMPS.  In  the  first  simulation,  only  the  atmospheric  model  was  used,  assuming 
a  steady-state  (fixed)  SST.  In  the  second  simulation,  both  the  atmospheric  and 
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ocean  models  were  used,  allowing  the  ocean,  and  therefore,  the  SST,  to  respond 
to  the  surface  fluxes.  In  both  simulations,  the  hurricane  remained  in  the  center 
of  the  grid  because  of  the  absence  of  a  mean  flow  and  use  of  an  f-plane.  The 
initial  thermodynamic  state  for  the  atmosphere  in  both  simulations  was  based  on 
an  individual  sounding  representative  of  the  environment  near  Gilbert.  A  Rankine 
vortex  was  used  to  initialize  the  hurricane  circulation.  A  climatological  temperature 
and  salinity  profile  was  used  to  initialize  the  ocean. 

For  the  fixed  SST  simulation,  the  model  generated  hurricane  showed  strong 
resemblance  to  many  features  observed  in  hurricanes.  The  central  pressure  fell  to 
874  mb,  with  a  deepening  rate  similar  to  that  found  in  Hurricane  Gilbert.  The 
model  storm  developed  an  eye  wall  with  a  radius  of  20-30  km  which  was  charac¬ 
terized  by  upward  motion  with  maximiim  values  of  up  to  6  m  s"^.  Subsidence  in 
the  eye  was  as  strong  as  -1  m  s"^.  After  an  initial  spin-up  period,  the  model  storm 
developed  spiral  bands,  consisting  of  convective  cells  which  propagated  in  toward 
the  eye  wall.  The  spiral  bands  would  occasionally  close  off  into  an  outer  eye  waJl 
which  would  contract  in  while  the  inner  eye  wall  weakened.  During  this  stage,  the 
storm  would  remain  at  the  same  intensity  or  weaken  slowly.  Eventually,  the  two  eye 
walls  would  merge,  and  in  one  instance,  this  wais  associated  with  a  central  pressure 
drop  of  nearly  10  mb  in  one  hour. 

In  the  coupled  simulation,  the  intensity  of  the  hurricane  was  hindered  by 
strong  cooling  of  the  SST  forced  by  the  hurricane  circulation.  Differences  between 
the  non-coupled  and  coupled  simulations  were  less  than  10  mb  for  the  first  42  h 
of  the  simulations  while  the  largest  SST  anomaly  was  no  more  than  0.5®C.  The 
lowest  pressure  reached  in  the  coupled  simulation  however,  was  914  mb,  which 
occurred  at  about  60  h  into  the  simulation,  or  about  12  h  after  the  SST  began 
to  cool  rapidly  near  the  center  of  the  storm.  After  60  h,  the  coupled  simulation 
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exhibited  a  weakening  of  approximately  15  mb.  After  84  h  of  integration,  a  SST 
anomaly  of  8°C  was  found  in  the  circular  area  bounded  by  the  radius  of  maximum 
winds.  The  central  pressure  of  the  non-coupled  and  the  coupled  cases  differed  by 
54  mb.  Some  enhanced  cooling,  of  up  to  0.75“C,  was  found  to  occur  under  the 
spiral  bands  that  formed  outside  the  eye  wall.  Under  these  bands,  the  near-surface 
air  temperature  was  found  to  be  lower  than  the  surrounding  air  due  to  evaporation 
of  raindrops.  This  led  to  enhanced  surface  sensible  heat  fluxes  and  stronger  SST 
cooling.  Finally,  it  was  found  that  although  the  SST  anomaly  of  8®C  altered  the 
intensity  of  the  system  near  the  eye  wall,  it  had  little  or  no  effect  on  the  outer  storm 
strength. 

The  following  conclusions  are  made  based  on  these  simulations.  First,  the 
feedback  of  ocean  cooling  on  the  hurricane  should  be  minimal  in  real  situations 
in  which  the  storm  has  a  translation  speed  of  at  least  1-2  m  s"*.  This  is  because 
of  the  approximate  12  hour  time-lag  of  response  of  the  atmosphere  to  the  SST 
changes,  and  that  during  this  12  hour  period  the  storm  will  i.iOve  to  an  area  with 
an  undisturbed  SST.  Second,  any  changes  that  are  made  to  the  tropical  cyclone 
due  to  an  induced  SST  coohng  will  be  most  evident  near  the  inner  eye  wall,  but 
not  in  the  outer  storm  strength.  This  would  imply  that  the  induced  SST  anomaly 
may  not  affect  the  track  of  the  system. 

The  results  presented  here  represent  only  a  fraction  of  the  type  of  studies  that 
will  be  possible  using  a  coupled  mesoscale  prediction  system.  First,  we  anticipate 
an  expansion  of  the  work  presented  here  to  include  moving  storms,  as  well  as  effects 
of  radiation  and  landfall.  Second,  we  are  working  toward  expanding  COAMPS  to 
include  topography  in  both  the  atmosphere  and  ocean,  radiation  parameterization 
and  modifications  to  allow  us  to  use  COAMPS  in  real-data  simulations.  Third, 
we  are  currently  validating  COAMPS  ability  to  simulate  other  features  such  as 
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rotating  convective  clouds,  topographically  induced  flows,  lake-eflect  snowstorms 
and  the  structure  of  the  planetary  boundary  layer  in  coastal  areas. 
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